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Abstract 

The effects of the intrinsic charm on the hadronic production of H cc are studied. By taking 
reasonable intrinsic charm component into account, the change of the theoretical prediction on 
the production of H cc for LHC and Tevatron is small, but in contrast it may enhance significantly 
for SELEX. The reason is that the collision energy at LHC and Tevatron is so large that the 
gluon-gluon fusion sub-process, which is irrelevant to intrinsic charm, becomes dominant. But the 
situation for SELEX is quite different. Our numerical results for SELEX show that by considering 
all the contributions from various sub-processes, the predicted cross-section may be enhanced by a 
factor so big as 10 2 due to a modulating intrinsic charm being taken into account. Therefore, the 
hadronic production of H cc at SELEX may be sensitive enough in observing the intrinsic charm 
inside the incident hadrons. 



PACS numbers: 12.38.Bx, 12.39.Jh, 13.87.Ce, 14.20.Lq 
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I. INTRODUCTION 



SELEX Collaboration [l, 2] has reported the observation of the doubly charmed baryon 
H+ which contains two charm quarks. According to their measurement, the decay width and 



production rate, nevertheless, are much larger than expected 



The theoretical 

predictions based on gluon-gluon fusions in [3|, |4], |5|, |6|, |7|] are much smaller than the measured 
cross-section by order of about 10 3 . Several possible mechanisms for the hadronic production 
of H+ are proposed, the discrepancy between experimental data and theoretical predications 

n 

decreases but it is still there [8( and it is at order of about 10 2 . Since the baryon S cc , either 
H+ or is doubly charmed, here we will highlight the intrinsic charm in the colliding 
hadrons through studying the hadronic production of S cc at different situations (e.g. at 
LHC, Tevatron and SELEX). 

In the hadronic production, the S cc is produced through scattering or annihilating or 
fusion of two initial partons inside the incident hadrons. Although the common belief is 
that the gluon-gluon fusion mechanism is important and even dominant in high energy 
hadron collision, the other mechanisms may be still substantial sometimes. For the other 
mechanisms, besides the light partons, it is also possible to find certain importance of heavy 
quarks inside the initial hadrons, whose possibilities are described by the corresponding 
heavy quark distribution function. The lightest heavy-flavored charm-quark, being as a 
parton, can be generated in different ways inside the incident hadrons. The heavy quark 
(charm or bottom) component in parton distribution functions (PDFs) can be perturbatively 
generated by gluon splitting, and hence is named as the 'extrinsic' component according to 
Ref.js]. It can also be generated non-perturbatively and appears at or even below the 
energy scale of heavy quark threshold. This component is normally called as 'intrinsic' one 
according to Refs. 9, 10, 11, Q, Q|- The extrinsic charm contribution to the S cc hadronic 



production has been studied in Ref.|8l within the general-mass variable-flavor-number (GM- 

nnn 

VFN) scheme \l4L Il5l . 1 1 61 ] . It is shown that the extrinsic charm contribution is substantially 



large at the energy range of SELEX experiment [8j, i.e., the extrinsic charm can raise the 
total hadronic cross-section of S cc by almost an order of magnitude in comparison with that 
only from the gluon-gluon fusion mechanism. It is also interesting to investigate how much 
the intrinsic charm may contribute in the hadronic production, specially to see where and 
to what 'degree' the intrinsic charm component in concerned hadron can be determined. In 
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fact, the existence of the intrinsic heavy flavor components in a proton has been supported 
by some experiments and theoretical predictions. For instance, the intrinsic heavy flavors 
are adopted to study the diffractive Higgs production in Ref.[17||. where the authors have 
pointed out that a clear experimental signal for Higgs can be observed due to the existence 
of intrinsic heavy flavor components inside the nucleon. Bearing the problem of the intrinsic 
charm inside a nucleon in mind, we will examine the intrinsic charm effects on the S cc 
hadronic production at different experiment situations, i.e. at LHC, Tevatron and SELEX, 
and see which is the best place to perform the measurements on the intrinsic charms 1 . 
In fact, the intrinsic charm quark content of proton is still not well determined. In 



the latest versions of PDF, like CTEQ [18j, GRV [19] and MRST [20], the heavy quark 
components arise only perturbatively through gluon splitting, which means they have the 
extrinsic nature. Although these PDFs are determined through the global fitting, it can not 
rule out the intrinsic component with certain possibility because the data for fitting does not 
include those experiment results which are sensitive to the intrinsic heavy quark reactions 
{^J. At present, the probability magnitude of the intrinsic charm within nucleon has not 
been well fixed yet. Nevertheless, we know it can not be be large in order to comply with 
the global fitting result. Upper-bounds for the probability of the intrinsic charm A in inside 
nucleon, i.e., the first moment of the intrinsic charm quark PDF at a lower energy scale (e.g. 



0,Q. 
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121 ]. A in can not be much larger than 



/i 2 ~ m 2 c ) is available: according to Refs 
1%; and another one, (0.86 ± 0.60)%, is obtained by analysis of data in deep inelastic 



muon scattering on iron 



21 



221 ]. The intrinsic charm PDF can only be studied with some 
nonperturbative methods with models. Recently, the operator product expansion method 
shows that the probability for Fock states in light hadron to have an extra heavy quark pair 



of mass Mq decreases as Aq CD /Mq in non-Abelian gauge theory 23j. In the present work 
for definiteness, we will take the BHPS model typical example for the intrinsic 

charm to study the S cc hadronic production. 

The paper is organized as follows. In Sec. II, we outline the techniques for calculating 
the S cc hadronic production under the GM-VFN scheme. Specially, we will show how we 
obtain the requested PDF of charm quark at the higher energy scale through CTEQ6HQ 



1 Since there are some puzzles in hadronic production of J/ip still, thus we think that any determination 
on the intrinsic charm components of the PDFs of the incident hadrons cannot be decisive. 
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plus the intrinsic charm component obtained by solving DGLAP equation with proper initial 
boundary condition at low energy scale. In Sec. Ill, we present the numerical results for the 
S cc hadronic production. The final section is reserved for discussions and summary. 



II. CALCULATION TECHNIQUES 

With QCD factorization and the GM-VFN scheme, the differential cross-section for the 
inclusive hadronic production of Hi + H 2 — > H cc + X can be written as the convolution 

+ { Yj fk (Xtoti-fHj (X 2 , 11)sUb\ ® Vgc->S cc (x 1 ,X 2 ,n)} 

+ { E [(/ffiC^l.M) - fuM^^SUB) ■ (f C Hj {x2,H) ~ f C Hj {x2,fi),SUB 



)}+•■■, (1) 



where \x is the renormalization (factorization) scale, the ellipsis stands for contributions from 
higher order in a s and from light quarks initiate processes. The later are much smaller than 
those given explicitly in Eq.(TJJ). ffj(a = g, c) is the PDF of the corresponding parton a and 
fff(x, fi)suB is defined in the GM-VFN scheme as 

r H {x^)suB = f g H (x,ri®ffa,n) = £jf c 9 (y,»)f 9 H (^») (2) 



with 



SiM - f^W^f^.^l-!,^. (3) 

To be realistic we will consider only the S cc with nonzero transverse momentum p t . 
Therefore, the partonic cross-sections in Eq.([T]) represent various parton processes at leading 
order in a s . The a gg -,~ cc stands for g + g — > E cc + c + c, cr gc ^ cc stands for g + c — > S cc + c 
and (Jcc-^Hcc for c + c — > E cc + g. These partonic cross- sections can be further factorized if we 
take charm quark as heavy quark. In this case, charm quarks inside the rest H cc move with 
a small velocity v. One can systematically expand the partonic cross-sections in Eq.(l) in v 
by using nonrelativistic QCD (NRQCD) {24] to separate perturbative and nonperturbative 



effects. In this framework the partonic cross-section can be expressed as 251 ]: 



a ab ^s cc = H(ab -> (cc^Sxjs) ■ h 3 + H(ab -> (cc)[ 1 5o] 6 ) ■ h + ■ ■ ■ , (4) 



where the ellipsis stands for the terms in higher orders of v. H(ab — > (cc) p\Si]g) or H{ab — > 
(ccj^Sole) is the perturbative coefficient for producing a cc pair in configuration of 3 Si 
and color 3, or 1 Sq and color 6 respectively |8j. The parameters /13 and /ii characterize 
the transitions of a (cc)[ 3 S , i]g pair and a (cc)[ 1 S'o]6 pair into the produced S cc , respectively. 
They are nonperturbative in nature and are just the relevant matrix elements in NRQCD 
framework (see Eq.(j3J)). It has been pointed out in 25| that these two parameters are at the 
same order of v, so for convenience and to decrease the 'freedom' of the prediction, hereafter 
we assume hi is equal to h 3 . 

Now we turn to the part on how to deal with the additional contributions from the 
intrinsic charm components to the charm PDF. It is generally expected that the intrinsic 
heavy quark component in PDFs is proportional to Aq CD /Mq with Mq being the mass of 
the heavy quark 23|. Therefore, in general the heavy quark components are small, so we will 
treat the intrinsic charm component as a small perturbation to the current determined PDFs. 
We will also take the leading order DGLAP equation for its and the relevant gluon component 
evolution. Hence, at leading order the introduction of the intrinsic charm component will 
not affect the PDFs of light quarks, but will affect the gluon's PDF. Under the above 
approximation the charm- and gluon- PDFs can be written as: 

/ff(^M) = fH^fav) + f C H n {xiV), (5) 
fij(x,fi) = ffi ex (x,fi) + f$ m (x,n). (6) 

In the above equations, fp ex (x,fi) is the extrinsic charm component that is already de- 
termined by the global fitting of several groups. fp m (x, /i) stands for the intrinsic charm 
component. The intrinsic charm component will also induce a small change to the currently 
determined gluon PDF, which is denoted as ffj m (x, fi). Hence, the gluon PDF is the sum 
of the change ffj m (x,[i) and the currently determined gluon PDF ffj ex {x, /i). 

The intrinsic charm component can not be calculated with pQCD. It can be introduced 
only with nonperturbative methods at some lower energy scale. Several models have been 
constructed for fp %n (x, //) M =2m c m a proton 0, The function fp* n (x, p) ll= 2m- from 

various models constructed in Ref . [{| is close to that of the BHPS model 13j] in shape, so we 
will take BHPS model as a typical one in our numerical studies. In this model the intrinsic 
charm component at \i — 2m c is parameterized as: 

fp in (x,2m c ) = fp in (x,2m c ) = 6£ [6a; (1 + x) lnx + (1 - x)(l + Wx + x 2 )] x 2 , (7) 
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where P stands for the proton, the parameter £ is determined by the first momentum of the 
distribution, i.e., the probability to find a charm quark in total: 



ft in (x,2m c ) dx = £x 1%. 



When £ = 1, it means that the probability for finding c/c-component in proton at the fixed 
low-energy scale 2m c is 1% as suggested in [3, Q]- In the following, we will take a broader 
range £ G [0.1, 1] to do our discussions 2 . The charm content in an anti-proton is the same 
as that in a proton. Since we will only deal proton-proton or proton-anti-proton scattering 
we will suppress the subscript H for PDFs and introduce shorthand notations for those in 
Eq.(4): 



c, in 
H 



in 

C 1 



rg, in 
J I 



fin rc, i 
Jq 1 JH 



_r0 r-g, ex _ rO 
J c ' J H in! 



f H ~ fci ffl ~ f< 



{I- 



(8) 



where H stands for a proton or anti-proton. 

With the intrinsic component fixed at the initial scale 2m c we can obtain the distributions 
at any higher scale n by employing the DGLAP equation. Taking the leading order equation, 
we use the the approximate method of Ref. [2 71] to solve the DGLAP equations for /* n (x, //) 
and fg n (x, /i). The solutions are: 



y { r(a/«) j 



dy f 1 dz 



i 

x y Jy z 
2 k 



f™(y/z,2m c 



\n(z)} 



an—l 



TidK) 



-Pac(x/v) + 0{k 2 ) 



'9 

with 



1 d Ji [" da f d -i /•«(,, 2m c 

Cb c J x y J a c Jy Z I 



8 2 (a s {2m c 
-, k = — in — — 

3 Po V a sW 



Hz)] 



CIK— 1 



T(aK) 



-P, 



c^gc 



x/y) + 0(k 2 1?) 



/5 = ll-2n / /3, 



Pac(x) 



P 



c^gc 



V^ + r- + (|- 2 7 B )^'-'-) 

1 — x mx V 2 



1 + (1 



x 



X 



(10) 



where n/ is the number of the flavor, and is taken to be 4. It should be noted that Pac is not 
exactly the splitting function P c ^ gc 27j. With the above equations we can obtain PDF's 



2 A quite smaller value of A. m ~ 10 5 has been suggested in Ref. 26|, since it is model dependent, we will 
not take such a small value to do our calculation. 
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FIG. 1: The ratio k c = f c (x, fJ>)/f® (x, fi) (left) and the ratio K g = f g (x,ij,)/fg(x,fi) (right). Two 
energy scales and three typical normalization for the intrinsic charm, Ai n = 0.1%, 0.3% and 1%, 
are taken respectively. 



at any scale. To make the employed approximation valid, it has been suggested in Ref. 



27] 



that k given in Eq.(9,10) should be smaller that 0.3. We have checked this in our case with 



/i = M t = \Jml cc + pi . For all ranges of \i used here k is smaller than 0.2. Since we include 
the intrinsic charm in PDF, the PDF of light quarks and gluons will also be changed in order 
to satisfy the momentum sum rule. This change can be expected to be at order of k 2 and 
leads to a violation of the momentum sum rule at 1% level with A in = 1%. We will simply 
neglect this change since the change is small. 

To see how the intrinsic charm component modifies the PDFs, the ratios k c = 
f c (x, (i)/ f®(x, fi) and K g = f g (x,fi)/fg(x,fi) are plotted in FigfTJ at the energy scales of 
fx 2 = lOGeV 2 and /i 2 = lOOGeV 2 with A in = 0.1%, 0.3% and 1%. The distributions with 
the same parameters are plotted in Figj2j From FigfTJ one can see that the intrinsic charm 
component changes the charm PDF significantly and k c increases as x increases. When 
x > 0.8 the change of the charm PDF after including the intrinsic charm component is too 
large to warrant the used approximations. Fortunately, the x which leads to the dominant 
contributions to differential cross section falls in the intermediate region. Our numerical 
study also shows that the PDFs in the region of x > 0.8 only give negligible contributions 
to differential cross-sections. From Fig.l and Fig.2 we can see that the change of the gluon 
PDF by including the intrinsic charm component is indistinctive as expected. 
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FIG. 2: The charm PDF (left) and gluon PDF (right) at fi 2 = 10 GeV 2 . The solid line is for 
xf®(x,fj,) or xfg(x,fi). The dashed line, the dotted line and the dash-dot line are for xf c (x,fj,) or 
xf g (x,fi) with Ai n = 0.1%, 0.3% and 1%, respectively. 

III. NUMERICAL RESULTS OF DIFFERENTIAL CROSS-SECTIONS 

In order to make numerical estimations on the differential cross-sections we need to know 
the parameters hi and h 3 . Unfortunately, they are unknown. For simplicity, we will take 
hi = h% as claimed in above. If one uses a nonrelativistic wave function ^/( cc ) for the 3 Si 



cc pair in the color 3 state, then we have h 3 = |^( cc )(0)| 2 [25|]. We take h% = |\l/( cc )(0)| 2 = 
0.039 GeV 3 as in [5[. The uncertainties from the values of hi and h 3 can be easily obtained, 
since they are overall factors for the processes. For masses we take m= cc = 3.50 GeV 
and m e J* = 1.75 GeV. The energy scale /i is fixed to be the transverse mass of H cc , i.e. 



fi = M t = \Jml cc + p 2 , where p t is the transverse momentum of the baryon. In dealing with 
the running coupling a s we take nf = 4 and Aq ( C^> 4 ' ) = 0.215 GeV. For all the numerical 
evaluations in the following we make use of the CTEQ6HQ version of PDFs. The analytical 



and numerical results are obtained by using well-established codes of BCVEGPY 
FDC 



and 



29j. Furthermore, one can convenient 



to study the hadornic production of S c 



y use the newly developed generator GENXICC 
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TABLE I: The contribution of a a b from different sub-processes initialized by the partons ab to the 
total cross section (in pb) for the H cc hadronic production at SELEX with the cut of pt > 0.2 GeV. 





CTEQ6HQ(Aj n = 0) 


A in = 1% 




°99 




°gc 


°99 






(cc)a[ 3 Si] 


4.03 


1.02 x 10" 3 


102. 


4.06 


1.25 x 10~ 2 


372. 


Wrf So] 


0.754 


4.15 x 10~ 5 


11.3 


0.758 


5.01 x 10" 4 


40.9 



TABLE II: The contribution rates of the sub-process gc — > H cc in the different x region in the 
charm quark PDFs with Ai n = 1% and pt > 0.2 GeV. 



0.0 <x c < 0.2 


0.2 <x c < 0.4 


0.4 <x c < 0.6 


0.6 < x c < 0.8 


0.8 < x c < 1.0 


25% 


50% 


22% 


3% 


~ 



A. Hadronic production of E cc at SELEX 

The cross-sections receive contributions from different sub-processes. To see how these 
sub-processes contribute to S cc production cross-sections we give our numerical results of 
different sub-processes in TABLEJTJ where we indicate explicitly the initial parton states 
and the final cc pair states. From TABLEjT] it is obvious that the intrinsic charm has the 
most significant impact on the contribution from the sub-process gc — > S cc . Also to note 
that this sub-process gives the dominant contribution to the total cross-section with A in = 
or A in = 1%. By taking A in = 1% the contribution is roughly four times larger than that 
without intrinsic charm. Taking the g + c sub-process as an explicit example, we present 
how the different regions of x in the charm quark PDFs contribute to the sub-process in 
TABLEjTTl From TABLEjTTl one can see that the region x > 0.6 gives tiny contribution 
and the main contribution comes from 0.2 < x < 0.4. We also note that in TABLEfl] the 
contribution from the quark pair cc in the configuration (ccJg^So] are generally smaller by 
a factor of about 10 _1 than that in the configuration (cc)g[ 1 S'o]. It should be noted that the 
contribution from (cc)g[ 1 ) So] is zero if we take hi = as indicated in Eq.(4). 

In Figj3]we plot the pt-distributions and the y(rapidity)-distributions, where contributions 
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FIG. 3: The ^-distributions (left) and y-distributions (right) for the hadroproduction of S cc at 
SELEX with different values of A^ n . The dotted, the dashed and the dash-dot lines are for Ai n = 
0.1%, 0.3% and 1% respectively. The result with CTEQ6HQ, i.e., Ai n = is shown by a solid line 
(the lowest one). 



TABLE III: The R values for SELEX with the cut p t > 0.2 GeV. 





CTEQ6HQ(yL in = 0) 


A in = 0.1% 


A in = 0.3% 


A in = 1% 


R 


29.3 


36.6 


51.3 


103. 



from all the above mentioned sub-processes are summed up. Our results show that the pt 
distributions do not change their shapes significantly, but the normalization changes. Also, 
as expected, the p t differential cross-section becomes larger as A in increases from to 1%. 
The shapes of ^-distributions and their normalization change significantly by changing A in 
from to 1%. With these results we may conclude that the intrinsic charm has significant 
impact on the production of S cc at SELEX. 

Before the work of Ref. [25| and without taking the charm PDF into account, the S cc 
production at a hadron collider is believed to be through the sub-process gg — > (cc^gpSi]. 
Now with the extrinsic/intrinsic charm and the (cc^^So] configuration contributions, there 
are additionally several other sub-processes which are non-negligible. To see how this alters 
the theoretical prediction based on the gg — > (cc)§[ 3 5i] subprocess, we introduce: 

R= ^ , (11) 

<7flfl-»Scc((cc) B pSi]) 

10 



2 4 6 8 10 12 14 

P t (GeV) 



FIG. 4: The energy scale dependence of the ^-distributions for g + c mechanism at SELEX. The 
upper band is for the case of (cc)jj [ 3 £i] and the lower band is for the case of (cc)e[ 1 5'o], where the 
solid line in each band corresponds to \i = Mf, the upper edge of the band is for fx = Mt/2 and 
the lower edge is for fi = 2Mt- Here the intrinsic charm has been taken into consideration, i.e. the 
PDFs of Eq.© with A in = 1% are used. 

where a to tai stands for the cross section with contributions from all sub-processes, including 
the intrinsic charm initiated ones; cr gg ->E cc ((cc)3[ 3 Si]) is the cross section only from the sub- 
process gg — > (cc)g[ 3 Si]. The numerical results of R for SELEX is given in TABLE JITII The 
results in TABLE JITT1 show that the cross section at SELEX will be enhanced by including 
all the dominant sub-processes contributions. The enhancement can be even up to orders of 
magnitudes, e.g. the prediction for the production cross-section can be enhanced by order of 
about 10 2 for the case of A in = 1%. Hence it is possible to reduce the discrepancy between 
theory and experiment from the order of 10 3 to 10 1 . 

Finally, we give some results of the uncertainty introduced by the factoriza- 
tion/renormalization scale. For clarity, we take the factorization scale Lip and the renormal- 



ization scale /j,r to be the same, /ijr = fiR = [i, and take three typical choices for /i 31] . i.e. 



LL — Mt (the default one in our calculations and M t = y M 2 + p 2 ), fj, = 2M t and fi = M t /2. 
The energy scale dependence of the pt-distributions for g + c mechanism are shown in FIGH] 
for SELEX. Numerically, one may find that by taking fi = M t /2 (or fi = 2M t ), the inte- 
grated cross-sections of the g + c — > S cc mechanism (either for the (cc)3[ 3 5i] configuration 
or the (cc) 6 [ 1 5'o] configuration) will be increased (or decreased) by about (2.0 ~ 4.0) times 
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10 15 

P t (GeV) 

FIG. 5: The pt-distributions for the hadroproduction of H cc at LHC. The left figure is for CMS 
or ATLAS with the rapidity cut \y\ < 1.5 being adopted and the right one is for LHCb with the 
pseudo-rapidity cut 1.8 < < 5.0 being adopted. The solid line, the dash-dot line and the circle 
line correspond to that of the g + g, g + c and c + c mechanisms without the intrinsic charm being 
considered (the PDFs in CTEQ6HQ [8] are used) respectively. The dotted line, the dashed line 
and the diamond line correspond to that of the g + g, g + c and c + c mechanisms with the intrinsic 
charm being considered (the PDFs of Eq.([6]) with Ai n = 1% are used) respectively. The differences 
with and without intrinsic charm are so small, that, of them, only at LHCb for the g + c mechanism 
the difference can be seen from the right figure. 

to the case of fi = M t within the allowable region of p t . 

B. Hadronic Production of E cc at Tevatron and LHC 

We have evaluated the S cc production at ATLAS or CMS and LHCb at LHC (the center 
of mass energy E cm = 14. TeV) and at CDF or DO at TEVATRON (the center of mass 
energy E cm = 1.96 TeV), where the intrinsic charm contribution is taken into account to 
see the possibility to observe the intrinsic charm effects at the high energy colliders. Since 
at LHC, the detectors ATLAS and CMS are similar, LHCb is quite different in rapidity 
cut; at TEVATRON the detectors CDF and DO are similar, thus we will present the results 
bearing the similarity and difference in mind of the detectors. We plot the ^-distributions 
in FIGs. ( 151161) with or without the intrinsic of Ai n = 1%, where the used cuts in rapidity 
are also given and the PDFs is taken from CTEQ6HQ [8] directly. From the figures one 



p (GeV) 
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FIG. 6: The pt-distributions for the hadroproduction of H cc at TEVATRON with the rapidity 
cut \y\ < 0.6 being adopted. The meaning for the lines in the figure is the same as FIGJ5J The 
differences between the two cases with and without intrinsic charm are too small to be seen. 





FIG. 7: The pseudo-rapidity 77-distributions for the hadroproduction of S cc through the g + c 
mechanism at LHC, where the pt cut pt > 4 GeV is adopted. The solid line, the dash-dot line 
correspond to the case of (cc)§ [ 3 5*i] and the case of (cc^^So] without the intrinsic charm being 
considered (the PDFs in CTEQ6HQ [8(] are used). The dashed line and the dotted line correspond 
to the case of (ce)g[ 3 5i] and the case of (cc^^So] with the intrinsic charm being considered (the 
PDFs of Eq.© with A in = 1% are used). 

can see that for all of the mechanisms g + g, g + c and c + c at the detectors of LHC and 
TEVATRON, the differences in ^-distributions for the two cases, i.e., with or without the 
intrinsic charm are too tiny to be seen in the most cases. The largest difference is found for 
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TABLE IV: Contributions from different sub processes to the cross sections for the hadronic pro- 
duction of S cc at Tevatron and LHC. The cut pt > 4 GeV is taken for all the hadronic colliders. 
Additionally, the cut \y\ < 1.5 is taken for LHC, the cut 1.8 < < 5.0 is taken at LHCb, and the 
cut |y| < 0.6 is taken at Tevatron. 





Tevatron 


LHC 


LHCb 








(cc) 5 [ 3 5i] 


Mel 1 So] 


(cc)3[ 3 Si] 




(Tgg (nb) 


1.61 


0.399 


22.3 


5.44 


25.7 


6.47 


(jg C (nb) 


2.31 


0.361 


22.1 


3.42 


20.6 


3.00 


a cc (nb) 


0.755 


0.0435 


8.75 


0.478 


3.18 


0.169 



TABLE V: The value of e as defined in Eq. (|12|) for the hadronic production of H cc at Tevatron and 
LHC. In the calculation, the cut \y\ < 1.5 is taken for LHC, the cut 1.8 < \rj\ < 5.0 is taken for 
LHCb, and the cut \y\ < 0.6 is taken for Tevatron. 





Tevatron 


LHC 


LHCb 


Ptcut 


10 GeV 


15 GeV 


20 GeV 


10 GeV 


15 GeV 


20 GeV 


10 GeV 


15 GeV 


20 GeV 


£gc (Ptcut) 


~ 1.1% 


~ 1.2% 


- 2.2% 


~ 0.1% 


~ 0.2% 


~ 0.5% 


~ 6.1% 


~ 10.% 


~ 15.% 


&cc (Ptcut) 


~ 1.5% 


~ 2.2% 


~ 4.3% 


~ 0.4% 


~ 0.5% 


~ 0.6% 


~ 0.8% 


~ 2.0% 


~ 3.5% 



the mechanism g + c at LHCb when p t is larger than 15GeV (see the right figure in FIG. ([5])). 
In order to see this more clearly, in Fig.fjTj), we plot the pseudo-rapidity distributions for the 
g+c mechanism at LHC with a p tcu t — 4 GeV. It can be found that the difference for the 
two cases is sizable only at large pseudo-rapidity region. 

In TABLE JlVl we show the contributions from different sub-processes to the cross sections. 
Our results are given with cuts for different colliders and are derived by taking Ai n = 1%. 
From our results one can see that at Tevatron and LHC both gg — > S cc subprocess and 
gc — ► H cc subprocess are dominant, unlike the situation at SELEX, where only the sub- 
process gc — > S cc is dominant. This implies that the intrinsic charm has no significant 
impact on the S cc production at the Tevatron and LHC. 
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FIG. 8: The pj-distributions for the hadroproduction of H cc . The dotted line, the dash-dot line, 
the circle line and the diamond line are those corresponding to LHCb, LHC, Tevatron and SELEX 
respectively with A{ n = 0, respectively The solid, the dashed line, the triangle line and the cross 
line are those corresponding to LHCb, LHC, Tevatron and SELEX with Ai n = 1%, respectively. 
Only at SELEX, the difference between the cases with and without intrinsic charm can be seen. 

To see how intrinsic charm leads to different contributions in different pt-regions, we 
introduce 



where i = cc and i = gc denote the contributions from the sub-processes cc — > S cc and 
gc — > H cc , respectively, of denotes the cross section without the intrinsic charm, while a 
denotes what with Ai n = 1%. The results for the ratio with different cuts are given in 
TABLE]V] From TABLE]Vl one can see that by introducing intrinsic charm, the maximal 
change in the cross sections under various cuts is only about 10%. This leads us to conclude 
that it may be impossible to observe the effects of intrinsic charm at Tevatron and LHC 
through the E cc production. 

Finally we give in FigJSJthe transverse distributions of S cc for different experiments. From 
FigJH] we see that only at SELEX the effect of the intrinsic charm is significant because the 
sub-process gc — > S cc is dominant one. At LHC and Tevatron, the sub-process gg — > H cc 
becomes more important, and can even be dominant with the increasing p t . This is the 
main reason why our results for LHC and Tevatron do not show any significant effect from 
the intrinsic charm. 




&i(pt > Ptcut) ~ Vjipt > Vtcut) 
(?i{pt > Ptcut) 



x 100%, 



(12) 
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IV. DISCUSSIONS AND SUMMARY 



We have studied the S cc hadronic productions in experiments at SELEX, Tevatron and 
LHC. In our analysis we include contributions from various sub-processes, especially con- 
tributions from the intrinsic charm. For the intrinsic charm component inside a nucleon, 
we adopt the BHPS model 13j]. Our results show that the intrinsic charm can have signif- 
icant impact on the S cc production at SELEX. The cross section with certain reasonable 
kinematic cuts can be four times larger than that without considering the intrinsic charm. 
In contrast, the effect of the intrinsic charm is small at LHC and Tevatron. The main rea- 
son is that at those colliders with large energy, the contribution from the sub-process of 
gluon-gluon fusion becomes bigger than that from the sub-process of gluon-charm scattering 
(see TABLE jlV|) . and the intrinsic charm only makes a small effect on the gluon PDF (see 
Figs fTirSl) . Therefore, in comparison with all of the environments considered here and if the 
intrinsic charm component in a nucleon is really restricted in the region A in = 0.1%-1% as 
for the BHPS model, then only at SELEX the S cc production is suitable for 'measuring' the 
intrinsic charm component inside the incident hadrons. This is one of our main results. 

In the literature the hadronic production of H cc is believed to be dominated by the 
production of a cc pair through gluon-gluon fusion and the cc pair is in a pSjJg state, and 
it has been found that the predicted cross section based on this sub-process is much smaller 
than the experimental observation of SELEX, the discrepancy is roughly about 10 3 [3j]. The 
transition probability of a cc pair in [ 1 5'o]6 configuration to S cc is of the same order in v as 
what the cc pair in pS'ijg configuration 25j, so that the extrinsic charm contribution can 
enhance the production cross section quite a lot 8j . Here we show that by taking all the 
contributions from the intrinsic charm into account, the prediction for the production cross- 
section can be enhanced by order of about 10 2 for the case of A in = 1%, the discrepancy 
hence is reduced. To obtain a final solution for the discrepancy a more detailed study is 
needed. 
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